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Abstract We report a first attempt to produce metastable hep solid helium below 
its melting pressure. A focused sound pulse is emitted along the oaxis of a mono- 
domain hep helium-4 crystal starting from a static pressure just above the melting 
pressure. The sound pulse is made as simple as possible with one negative and one 
positive swing only. Density at focus is monitored by an optical interferometric 
method. Performed numerical simulations show that the crystal anisotropy splits 
the focused wave into two separate pulses, corresponding to a longitudinal wave 
along the oaxis and a radial one perpendicular to it. The amplification factor due 
to focusing remains nevertheless important. Negative pressure swings up to 0.9 
bar have been produced, crossing the static melting pressure limit. Improvements 
in the detection method and in the focusing amplification are proposed. 
PACS numbers: 67.80.B-, 67.80.bd, 64.70D-, 62.30.+d 



1 Introduction 

The hep solid phase for helium-4 is stable only above a minimum pressure Pf(T)~ 
25 bar at low temperature. Recently it has been suggested 1 that this solid phase 
remains metastable at much lower pressures. Based on an extrapolation of the 
equation of state, the spinodal line is predicted to be at negative pressure, about 
-10 bar. Beside testing this prediction, there is another interest to investigate solid 
helium at pressures below the melting pressure. Vacancies in this quantum solid 
have been studied for a long time. In particular their energy E v has been measured 
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at different pressures 2 , and also computed by Quantum Monte Carlo simulations 3 . 
Although there is a significant scatter on the measured values of E v , it is striking 
that it is a strongly decreasing function of the molar volume V m . A simple linear 
extrapolation on the data gathered in reference 2 suggests that an increase of V m by 
only 5% beyond the melting volume could bring E v near zero. Although this is a 
large increase for a condensed phase, this is not so for solid helium which is a very 
compressible material. More precisely, the maximum molar volume for stable hep 
helium-4^ is Vy = 21 cm 3 , the aimed value is V a = 22 cm 3 , while the predicted 
value of V m at the spinodal line is V s = 34.7 cm 3 . Hence V a is within a possibly 
accessible range. If so, this would bring interesting new physics for solid helium. 
With E v near zero, vacancies could proliferate, reach a sizable population, possi- 
bly undergo a Bose-Einstein condensation, thus realizing the Andreev scenario for 
supersolidity 5 . Or if the vacancy-vacancy interaction is strong enough— it could 
provide another mechanism for destabilizing the solid phase before reaching V s . 

Many different methods have been used to study metastable condensed phases, 
mostly liquid, under depression 7 . Metastability cannot be easily obtained for a 
solid because the interfaces between a solid and the container walls are disordered, 
and generally nucleates the liquid phase. Focused sound waves provide a way to 
circumvent this difficulty. Starting from a pressure Pq above Pt, a sound pulse 
with an amplitude 8P, smaller than Pq — Pf at the surface of the transducer will 
not bring the solid below the melting pressure at the transducer surface. On its 
propagation to the focus the sound pressure can be amplified by a large amount, 
namely Q.R/X S (Q, solid angle of the sound beam, R transducer radius, A s sound 
wavelength). Hence the peak pressure Po + 8PjQR/X s can explore the metastable 
domain below Pt with negative 8 Pi. 

This article reports a first attempt to produce metastable hep solid helium be- 
low its melting pressure. A hemispherical transducer was used to produce a con- 
verging sound wave. The density variations at focus are monitored by an optical 
interferometric method. The experimental arrangement is described in section 2. 
Because the sound velocity is anisotropic, the focusing is not expected to be per- 
fect in hep helium. A numerical simulation was made to model the remaining 
amplification in this case and to determine the relation between the optical signal 
and the density at focus. Numerical predictions are discussed in section 3. Experi- 
mental results are reported and discussed in section 4. Possible improvements and 
conclusions are given in section 5. 

2 Experimental arrangement and procedure 

The experimental cell is a 4 cm stainless steel cube, with five silica windows 
(diam. 2.5 cm), cooled from the top by a pumped helium-4 fridge in the 1.0- 
1.4 K temperature range. To achieve single crystal growth, the nucleation and 
growth is made at constant temperature and pressure (1.2 K and 25.5 bar). An 
electro-crystallization device 8 provides a unique seed which falls on the bottom 
of the cell. The crystal is subsequently grown from this seed. The orientation of 
the c-facet (most often horizontal) is easily monitored visually during the growth 
process which takes place below the corresponding roughening transition 1.3 K— . 

A hemispherical piezoelectric transducer of 6 mm internal radius (from Chan- 
nel Industries) is suspended above the middle of the cell with its axis vertical, 
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along the crystal oaxis. Such a transducer has two main resonant modes: the 
thickness mode frequency is V/,/ ~ 1 MHz, the inner and outer surfaces vibrat- 
ing in opposite directions. The breathing mode corresponds to an oscillation of 
the transducer radius and has a frequency V//- ~ 180kHz. In order to simplify as 
much as possible the pressure wave, a single oscillation of the thickness mode 
was used. The driving voltage is produce by an arbitrary function generator and 
amplified by a RF amplifier up to a voltage ranging from 100 V to 800 V on 50 Q. 
The pulse shape is designed to leave the transducer at rest at the end of the pulse. 
Beside the one-cycle oscillation at 1 MHz, the exciting pulse also produces a par- 
asitic oscillation of the low frequency breathing mode which appears after the end 
of the 1 MHz pulse. The sound velocity in hep solid helium is anisotropic. For 
longitudinal waves, the velocity is about 540 m/s along the c-axis and 460 m/s 
when the wave vector is in the basal plane 10 . Hence the initially spherical wave 
does not remain so and breaks up into two parts: one along the c-axis, named as 'z- 
pulse' hereafter, and another radial one, named V-pulse'. Both show a maximum 
amplitude near the center of the transducer. Their relative intensity and shape will 
be discussed in the next section. In order to have an optical access to the transducer 
center two small notches (1 mm wide, 1.5 mm high) were made on the transducer 
rim, along a diameter. 

The phenomena at the focus are monitored optically. A CW-laser beam (A = 
532 nm) propagating along the y-axis is focused at the cell center with a 30 fim 
waist. The laser polarization is along the vertical z-axis. The axis origin is taken at 
the transducer center. Density modulations along the beam path produced by the 
sound wave result in an optical phase change 



where / is the cell length. The change in refractive index 8n is related to that 
of density 8p through 8n/(n— 1) = 8p /p. This linear approximation of the 
Clausius-Mossotti relation holds for helium with a relative error less than 10~ 4 
for 8p/p up to 0.1, well beyond the strains discussed here (less than 10~ 2 ). Be- 
cause of the system cylindrical symmetry around the z-axis, 8n and 8p only de- 
pend on r = x 2 +y 2 and z. To measure <5</>(f), the laser beam is split before 
the cell to produce a reference beam which crosses the cell in an unperturbed re- 
gion and is recombined with the monitoring beam on a photo-detector, making 
a Jamin interferometer. 8(j)(t) contains a low frequency part 8(j>if(t) due to the 
transducer breathing mode (at V/f), and the interesting signal <5 </>/,/■(?) due to the 
thickness oscillation (at V/,/). To extract the latter from <5</>(f), a model of the low 
frequency signal As/f(t — to) has been constructed with two adjustable parameters: 
its amplitude A and its starting time to- These parameters are determined from a 
fit of As[f(t — to) to 8(j)(t) filtered with a low-pass filter. Then 5</>/,/-(?) is found as 
8<j) h f(t) = 8<l)(t)-As lf (t-to). 

Extracting 8n(t) from <5</>/,/-(f) is more involved. This amounts to determine a 
characteristic length a of the focused sound wave such that 




(1) 



5^/(^0,0) =2nSn(t,0,0)a/X, 



(2) 



For a continuous sinusoidal wave, a = X s j7} 2 . There are two ways to get a. If one 
can measure 8(j>hf(t,x,0) for the whole range of x where it is non-zero, one can 
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invert formula[T]and get Sn(t,x,0) from S(j>hf(t,x,0) by an inverse Abel transform. 
This is not always possible however due to the small width of the notch. A simpler 
method relies on numerical simulation to get the ratio <5</>/,/-(?,0,0)/n(?,0,0) and 
thus a. In the following, results are discussed using the optical length variations 
8L opt = X 8tyhf / (2n) as the signal instead of 8(phf- 



3 Numerical simulation of a focused spherical sound wave 

Since solid helium has a close-packed hexagonal structure (hep), its linear elas- 
tic properties are described by five elastic constants, c\\, cyi, c^, C33, C44 as 
defined by Musgrave 13 . The ratio of shear stress to shear strain in the meridian 
plane, called C55 is not an independent constant in hep crystals 14 and is equal to 
(ci 1 —C]2)/2. Numerical values of the c (/ have been taken from™: c\\j p = 21 .2, 
cn/p = 11.1, C13/P = 5.49, c 33 /p = 29.0, and C44/P = 6.52, all in units of 
10 4 (m/s) 2 . 

The axis of the transducer being parallel to the oaxis of the crystal, deforma- 
tions are of cylindrical symmetry. Elastic waves are thus governed by two coupled 
second order differential equations for the only two possible components of the 
displacement vector, u r and u z : 

d 2 u r \du r u r d 2 u, d 2 u r 

pii r = c\\ (-5-5- H 5 2) + ( c l3+ c 55)3-3-+c 5 5 3-^- (3) 

dr A r or r L drdz dz 

, ..d 2 u r 1 du r d 2 u, ,d 2 u z 1 du, . 

pu z = (ci3+c 55 )(3-3-H 5-)+ c 33 3-7-+c 5 5(3-^-H 5-) (4) 

oraz r dz dz dr L r dr 

(C44 plays no role in this geometry) 

Special care must be taken in r = 0: for symmetry reasons, u r , u r and u r must 
be in r = 0. Using the 1' Hospital's rule, equation|4]can be rewritten: 

„. . d 2 u r d 2 u, d 2 u z , 

pu z = 2(ci 3 +c 55 )-^-+c 33 -^ 1 -+2c 55 -^- (r = 0) (5) 

These equations are integrated using a finite difference method and a stag- 
gered leap frog scheme with initial conditions corresponding to an undeformed 
crystal at rest. The model volume is < r < 7.5 and —7.5 < z < 7.5 (in mm) 
and is discretized on a 1024 x 2048 grid. The spatial step 81 is then small enough 
compared to the wavelength (A s ~ 0.5 mm) that the dispersion introduced by the 
discretization is negligible. The time step 8t is taken to be sufficiently small to ful- 
fill the Courant criterion: essentially one must have 81/ 8t < X s Vf,f. At each time 
step t one then calculates a new field of displacement vectors with the help of the 
above mentioned differential equations as a function of the fields at time t — 8t 
and t — 28t. The field is then adjusted so as to fulfill the boundary conditions: u r 
and u z are forced to be that of the transducer on the points belonging to it and 
u z (r — 0,z) is adjusted to the value u z (0,z) — (4u z (8l,z) — u z (28l,z))/3 to ensure 
cancelation of its derivative (r = 0, z) ■ The other boundaries are left free, giving 
rise to reflected waves, which arrive at the focus with a long enough delay to be of 
no concern. From the displacement vector field, the change in molar volume can 
be computed (8V m /V m = + y + and thus the change in refractive index 
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Fig. 1 Map in the zx-plane of the refractive index variations while the mono-oscillation sound 
wave is propagating. The thick black line represents the section of the transducer surface. The 
initially hemispherical wave splits into an axial wave and a radial one. 



8n. An example of the wave resulting of a single oscillation of the transducer is 
shown on FigureQ]at time t = 9 /is. As expected, the wave has split into two parts, 
one along the z-axis and a slower one which becomes almost cylindrical, traveling 
roughly in the r direction. At the center one thus expects two pulses, as confirmed 
by Figure [2}a. When passing through the center, the r-wave undergoes a Gouy 
phase reversal. As a result the r-pulse is tripolar instead of being bipolar like the 
z-pulse. Note that this is true only at the center, because only at the center does 
the incoming and the outgoing waves interfere. An important result is the ampli- 
fication factor. The z-pulse is almost not amplified, while the r-pulse is amplified 
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Fig. 2 (a) Refractive index modulation near the transducer surface (dashed line) and at the center 
(solid line). The small modulation around 11 /is is the part of the sound wave coming along the 
z-axis (z-pulse). The large one around 13.5 /is is the radial one (r-pulse). (b) Modulation of the 
optical length through the center of the transducer for the same pulse. 
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by a factor 50. This is not much less than the factor 75 expected for an isotropic 
medium with the same X s . Figure [2] also shows the behavior of the optical length 
through the center of the transducer. A first striking result is that the optical length 
is vanishingly small for almost all times, although the integration path crosses the 
wave at any time. This is because the positive and the negative parts of the wave 
cancel almost exactly. Only when the wave reaches the center is this no more true. 
The z-pulse seems to be magnified, but this is only a geometrical effect due to 
the fact that this part of the wave is almost plane and perpendicular to the z-axis. 
The ratio between the maximum amplitude of the optical length and the maximum 
amplitude of the index modulation is a characteristic length a {a ~ 0.28 mm). It 
can be used to deduce the refractive index modulation, which is the quantity of 
interest, from the optical length, which is the measured quantity. 



4 Presentation and discussion of some results 

An example of observed optical signal <5L opt (/) is shown in Figure [3}a. It was 
taken with a driving voltage Vd = 200 Vatr=1.16K and a pressure 0.6 bar above 
Pf{T). An order of magnitude of the sound amplitude at the transducer surface can 
evaluated from the formula dp/p = 2Jtdy}Vd / X s f r , where = 0.3 nm/V is the 
piezoelectric constant at room temperature of the transducer ceramic, and f r ~ 4 
is an order of magnitude of its reduction factor at low temperatures 16 . One finds 
8p/p ~2x 10" 4 , which can be converted to an acoustic pressure dP = 54 mbar 
using solid helium bulk modulus B s = 269 bar. 17 . Thus the minimum pressure at 
the transducer surface is well above the melting pressure. Note that for low sound 
amplitudes (Sp jp < 10~ 2 ), the computation of 8P using the bulk modulus instead 
of a non linear equation of state does not lead to a relative error larger than 10~ 2 . 
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(a) (b) 

Fig. 3 Observed optical length modulation by focused sound pulses. Time origin is the begin- 
ning of the transducer oscillation which lasts 1 /xs. Arrows indicate the middle of expected z- and 
r-pulses. (a) driving voltage 200 V, initial offset from melting pressure 0.6 bar, pressure swing 
at focus -0.35 bar. (b) Same signal with parameter values 600 V, 0. 1 bar, -0.9 bar respectively. 

Comparison of the optical signal with Figure |2]-b helps to identify the radial 
wave. The expected time of flight is indicated by the arrow labeled 'r-pulse'. Note 
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that in Figures [2] and [3] the origin of time is the beginning of the transducer oscil- 
lation. To get the time of flight for the middle of the sound pulse, half its duration 
is to be subtracted, namely 0.5/is. The sound velocity at this particular pressure is 
determined using the scaling of c\\ and C33 proposed by Maris 1 . For this signal, 
the laser beam is thought to cross the focal region. In fact extensive scanning of 
the focal region along the x and z-axis with the laser beam was performed to locate 
the sound focus. 

The shape of the /--pulse is qualitatively similar to the expected one. The width 
in time measured by the delay between the secondary maximums is close to the 
1.0 /is delay in Figure Ub. There are several differences however. Some are par- 
asitic effects. The high frequency modulation may be attributed to higher oscilla- 
tion modes of the transducer which are found near 2.4 MHz and 3.6MHz. Oscil- 
lations before the expected z-pulse can be explained by the fact that the reference 
light beam is slightly perturbed by the wave emitted by the rim of the piezo- 
hemisphere.The distance between this rim and the reference beam is 4 mm and 
thus the sound wave reaches it in 7 /is. Concerning the /--pulse itself, the ratio of 
the main peak amplitude to the secondary maximums is much smaller than in the 
simulation. This could be due to an offset of the laser beam from the radial geom- 
etry or spatial averaging over the laser waist. Simulation shows that an offset of 
0.1 mm would divide the peak by a factor 3. 

There is not a clear evidence for the z-pulse. It is either absent or hidden by 
parasitic oscillations. Its weakening may be correlated to that of the /--pulse peak. 
Another possible explanation is that the initially well oriented mono-crystal has 
been broken and turned to a more or less isotropic polycrystal by the sound wave. 
As a matter of fact, during long measurement series for scanning the laser beam 
position, we occasionally observed evolution of the optical signal, generally a 
weakening of the pulse amplitude. 

Let us come back to the intensity of the negative swing at focus. Using the 
value of a determined by the simulation (a = 0.28 mm), one gets an order of 
magnitude of the negative pressure swing from the formula dP = 8L opt B s /a(n — 
1 ) ~ —0.35 bar. The initial pressure was 0.6 bar from the melting pressure, so that 
the minimum pressure does not reach the melting line. Figure [3}b shows another 
signal taken with a larger driving voltage (600 V) and starting closer to the melting 
line, P ~Pf = 0. 1 1 bar. The signal is clearly more distorted than in Figure [3]-a. If 
one converts nevertheless the negative swing into pressure using the same formula 
quoted above, one gets 8P = —0.9 bar. In that case, it appears that the pressure 
has crossed the equilibrium melting line by at least 0.8 bar. 

Has helium remained in the solid phase? Had a liquid droplet nucleated at the 
focus, it would have expanded for a fraction of the sound half period at a velocity 
close to the sound velocity. Thus it would have reached a diameter on the order of 
100 fim. Taking into account the refractive index difference between liquid and 
solid phases (about 0.0035), the liquid bubble would have created a jump in 8L opt 
as large as 0.35 /im. Also a large asymmetry would be found between positive and 
negative pressure swings. Neither of the two phenomena are observed and it can 
be concluded that transient melting did not occur. 
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5 Possible improvements and conclusions 

A direct interferometric imaging has been developed 18 to provide a map of 8L opt . 
It uses the same interferometer, but the laser is a pulsed laser and the detector a 
CCD camera. Images are taken for various / and reference phases. Then a map 
of 8(j>(t,x,z) can be computed for all the pixels in parallel. Beside providing a 
direct image of the sound pulse propagation in the crystal, inverse Abel transform 
will give access to the map of 8n(t,r,z), provided that it has a cylindrical sym- 
metry. This should overcome the major shortcoming in the present experiment, 
namely the conversion from the measured 8L opt to the pressure swing 8P at the 
focus. Since the observed signal differs from the simulation, a more direct way 
to computed the latter will be valuable. It will also avoid scanning the laser beam 
to locate the focus, which is both time consuming and harmful for the crystal. In 
order to view the focus and an extended portion of the x-axis, the lower part of 
the hemispherical transducer will be shorten by 1 mm or so. This new method will 
also help assessing whether the crystal is broken by repetitive sound pulses. If this 
is the case, it would be of interest to try to work with a polycrystal with grains 
smaller than k s . The medium will then be quasi-isotropic 19 . Whether the sound 
attenuation will be acceptable remains an open question. In the case monocrystals 
survive their repetitive stress, it would be more efficient to use a transducer with 
a shape conforming the wavesurface of the hep crystal, or at least half of it. As 
depicted in reference^, it is an elongated bowl, nearly conical in the directions 
between z and x-axis. Finally using longer multi-oscillation excitation pulses will 
allow to reach oscillation amplitudes 8P about 10 times larger. The more complex 
wave pattern will not be a problem if inverse Abel transform can be performed. 

To summarize, it has been shown that interesting motivations exist to inves- 
tigate solid helium below the melting pressure. Numerical simulations confirmed 
that sound waves can be efficiently focused in the anisotropic hep crystal and pro- 
vides a tool to investigate the metastable region below the melting curve. A first 
example of incursion in this domain has been reported, and perspectives exist to 
enlarge the explored domain. 
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